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its scaling and regularity properties. We show that such a representation can be 
used to accurately determine it from the knowledge of the few first coefficients 
of the expansion for small magnetization. Revised values of amplitude ratios are 
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method from the e = A — d expansion. 
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1 Introduction 

Renormalization group (RG) theory of second order phase transitions provides 
a complete description of aU universal quantities (for a general reference on this 
article see for instance |jl|]). Among them the most studied are critical exponents, 
because they are easier to calculate, and because they have been used to test 
RG predictions by comparing them with other results (experiments, high or low 
temperature series expansion, Monte-Carlo simulations). For the 0{N) vector 
model, calculations are based upon the {(j)^)^ field theory 

^(</') = / {l[dM^)f + l>^2(/>\x) + ^X4c/>\x)f} d'^x . (1.1) 

We recall that near the critical temperature Tc X2 is a linear measure of the 
temperature. If we denote by A2C the value for which the theory becomes massless 
(T = Tc) then the parameter t 

t = X2-X2c(xT-Tc, (1.2) 

characterizes the deviation from the critical temperature. 

To determine exponents two strategies have been used. One follows Parisi's 
suggestion [0 and is based on perturbation series calculated directly in three 
dimensions. Series up to six loops have been generated 0] which have been 
summed by summation methods based on a Borel transformation (first estimates 
for = 1 were reported in [Q). The more accurate estimates are obtained after 
a conformal mapping which takes full advantage of the large order behaviour 
analysis For a review see |jl|,0. In this article, as a test of the variant of 

the summation method we use, we have reexamined their determination. 

Alternatively the e = 4 — d expansion known up to five loop order ^ has 
been summed by similar techniques 0] (note however that the authors used the 
slightly erroneous fifth order of Gorishny et at). Both expansions lead to con- 
sistent results. Moreover RG values for critical exponents have now resisted for 
many years to confrontation with lattice calculations and experimental determi- 
nations. Table 1 displays results for Ising-like systems (A^ — 1) coming from field 
theory under the form of summed d = ?> series and e-expansion compared to ex- 
periments: binary fiuids (a), liquid-vapour transitions (6) and antiferromagnets 
(c), see [0,0 for references. Only recent Helium superfiuid transition experi- 
ments (A^ = 2) [|ll| now provide results which are consistent but more accurate 
than RG estimates. It would therefore be interesting to examine whether the 
series could now be extended to further decrease the apparent errors in d = 3 
dimensions. 

Other available estimates come from the analysis of High Temperature series 
in lattice models (table 2) and Monte-Carlo simulations. For the latter case let 
us quote two typical set of results in [l^] the values p — 0.629 ± 0.003 and r] — 



3 



Table 1 

Critical exponents for Ising-like systems: RG and experiments. 





7 


V 


a 


Q 


Q = ujv 


d = 3 


1.2405 ±0.0015 


0.6300 ±0.0015 


0.325 ±0.0015 


0.110 ±0.0045 


0.50 ±0.02 


£ — exp. 


1.2390 ±0.0025 


0.6310 ±0.0015 


0.327 ±0.0015 


0.110 ±0.0045 


0.51 ±0.03 


(a) 


1.236 ±0.008 


0.625 ±0.010 


0.325 ±0.005 


0.112 ±0.005 


0.50 ±0.03 




1.23-1.25 


0.625 ±0.006 


0.316-0.327 


0.107 ±0.006 


0.50 ±0.03 


(c) 


1.25 ±0.01 


0.64 ±0.01 


0.328 ±0.009 


0.112 ±0.007 





Table 2 

Critical exponents for Ising-like systems: HT series. 





7 


V 


a 


Q = uju 


1 
1 
1 

1 
1 
1 
1 


12 
13 
14 
15 
16 
17 
18 


1 
1 
1 

1 
1 
1 
1 


1.239 ±0.002 
1.2385 ±0.0025 

1.239 ±0.003 
1.2395 ± 0.0004 

1.239 ±0.003 

1.237 ±0.002 


0.631 ±0.003 
0.6305 ±0.0015 
0.631 ±0.004 
0.632 ±0.001 
0.632 ±0.003 
0.630 ±0.0015 


0.105 ±0.007 
0.101 ± .004 

0.104 ±0.004 


0.57 ±0.07 
0.54 ±0.05 
0.52 ±0.03 



0.027±0.005 are reported. In one finds u 
while 9 varies in the range 0.44-0.53. 



0.625±0.001 and ry = .025±0.006 



Other universal quantities, like the equation of state [pi| , |22| , |23 



and some am- 



plitude ratios 



have also been calculated but the estimates are less 
accurate because the series are shorter. Perturbative calculations are more dif- 
ficult, in particular for N ^ 1 (continuous symmetries), due to the presence of 
two different masses (transverse and longitudinal) in an external magnetic field 
and Goldstone singularities on the coexistence curve {H = 0, T < Tc). Probably 
also less effort have been invested up to now. 

In this article we therefore present a determination of the equation of state for 
the N = 1, d = 3 case, where longer series are available. Our calculations are 
based on perturbative expansion at fixed d = 3 dimension . Five loop series for 
the renormalized effective potential of the (p"^ theory have been first reported by 
Bagnuls et al. |2^, but the printed tables contain some serious misprints. These 



have been noticed by Halfkann and Dohm who have published corrected values 
(for detailed explanations see our appendix). Finally full analytic results up 
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to three loops have recently become available [p8 



The main technical difficulty that one faces m d = 3 calculations is how to 
continuously extrapolate field theory results from T > Tc to T < Tc. Indeed 
because the massless (or critical) theory is IR divergent in perturbation theory 
at any fixed dimension d < 4, calculations can be performed only in the massive 
phase (in contrast with the e-expansion) . In one method is suggested which 



has then been used [0 (see also [|6l) to predict some amplitude ratios. We 



present in this article a different approach that extends a suggestion in , and is 
motivated by the simplicity of the parametric representation |3^] of the equation 
of state within the framework of the e-expansion [^|23|| . 

The first step of our approach is a summation of the available perturbative 
series to obtain a non-perturbative determination of small field expansion of the 
effective potential (coefficients of (j)^, (p^, (j)^^) of the continuous cj)^ field theory at 
the IR fixed point. This result is interesting in itself in view of the recent ef- 
fort devoted to the problem (Monte Carlo lattice simulations, High Temperature 
series and approximate numerical solutions of Exact Renormalization Group). 
With the help of a parametric representation we are then able to reconstruct the 
full scaling equation of state, which is therefore the main result of our article. 
We then use it to calculate several amplitude ratios which have been considered 
before in the literature. We compare our predictions with other available results 
(e-expansion. High Temperature Expansion, Monte Carlo, Exact Renormaliza- 
tion Group and of course experiments). 

The set-up of our article is thus the following. In section 2 we briefiy review 
the properties of the equation of state and the definitions of amplitude ratios. 
In section 3 we recall known results about the e-expansion, while in section 4 
the idea of parametric equation of state and its application to the e-expansion 
are presented. In section 5 we explain the method we have used in our d = 3 
calculations. In section 6 the summation of the perturbative expansions for the 
effective potential is discussed. The results concerning the effective potential 
and the (parametric) equation of state are reported in section 7, while section 8 
contains our results for amplitude ratios and some concluding remarks. 



2 Effective action and equation of state 



In this article the general framework is the massive theory renormalized at zero 
momentum. The correlation functions T^^^ of the renormalized field (j)^- = ^l\fZ 
are fixed by the normalization conditions 

(p; m, £?) = + p2 ^ O (/) , (2.1a) 
Vf\v^ = ^■.m,g) = m^-''g. (2.16) 

In this renormalization scheme one trades the initial parameters A2 and A4 for g 
and m. The renormalized coupling constant g is dimensionless. It has to be set to 
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its IR fixed point value g* . In the value g* = 23. 73 ±0.09 has been numerically 
estimated from the series published in . The mass parameter m is proportional 
to the physical mass, or inverse correlation length, of the high temperature phase. 
It behaves for t oc T ~ Tc 0+ as m oc , where u is the correlation length 
exponent (see []T| for details). Therefore the parameter m is singular at and 
thus the extrapolation to the low temperature phase is non-trivial. Within the 
framework of the ^-expansion instead, it is possible to first construct the massless 
theory and then the theory in the full critical domain as an expansion in powers 
of the deviation t from the critical temperature. Equivalently the problem can 
be solved because scaling laws are exactly satisfied. In the summed d — 3 
perturbative expansion, at the IR fixed point g*, scaling relations are satisfied 
only within summation errors. The parametric representation of the equation of 
state will provide a solution to this problem. 



From the conditions it follows that the free energy J-" per unit volume 

expressed in terms of the "renormalized" magnetization 99, i.e. the expectation 
value of the renormalized field (py = (p/ -\/Z, has a small (f expansion of the form 
(in d dimensions) 

^(VP) = r((^)/vol. = ^(0) + imV' + im^-'^g^'' + O (cp^) , (2.2) 

where T{ifi) is also the generating functional of One Particle Irreducible correla- 
tors restricted to constant fields, or the effective potential of the renormalized 
theory. 

The equation of state is the relation between the magnetic field H, the mag- 
netization M = {(f)) (the "bare" field expectation value) and the temperature 
which is represented by the parameter t (eq. (|1.2|)) 

dM ' 

Near the critical point the equation of state has Widom's scaling form 

if(M, t) = M^f{t/M^'^). (2.3) 
It is convenient to introduce the rescaled variable z 

<^ = ^(^-2)/2^/^, (2.4) 

and set 

T{ip)-J^{Q) = —V{z,g). (2.5) 
9 

The critical behaviour is described by the IR fixed point, but we keep here the 
notation g to remind that we must first sum up perturbative expansion for V 
and then take g = g*. 
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The equation of state is obtained from the derivative F of the reduced effective 
potential V with respect to z 

= (2.6) 
Ising-like symmetry implies that -F has an expansion of the form 

F{z, g)=z + \z^ + Y. F2i+i{9)z''+' . (2.7) 



From eqs. ( |2.4| , |2.3|) we conclude 

z oc Mf-f^. (2.8) 

Comparing then the coefficient of M in the small M expansion of H{M, t) (see 
also eq. ( |2.13|) ) 

H{M,t) = Hof^M + (M^) , (2.9) 

where Hq is a numerical constant, with the small z expansion of the function 
F{z) we conclude 

H oc ff^^Fiz), (2.10) 
where the relation between exponents 

7 = /3((5-l), (2.11) 

has been used. One property of the function H{M, t) which plays an essential 
role in our analysis is Griffith's analyticity: it is regular at t = for M > fixed, 
and simultaneously it is regular at M = for t > fixed. 

Amplitude ratios. Universal amplitude ratios are numbers characterizing the 
behaviour of thermodynamical quantities near Tc in the critical domain, which 
in addition do not depend on the normalizations of magnetic field, magnetization 
and temperature. Several amplitude ratios commonly considered in the literature 
can be derived from the scaling equation of state: 

The specific heat. The singular part of the specific heat, i.e. the 4>'^ 2-point 
correlation function at zero momentum, behaves like 

CH = A^\t\~'^, tocT-Tc^±0. (2.12) 

The ratio A'^ /A~ then is universal. 

The magnetic susceptibility. The magnetic susceptibility x in zero field, i.e. the 
(p 2-point function at zero momentum, diverges like 

X = C^\t\~^, t^±0. (2.13) 
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The ratio /C~ then is also universal. 

Other amplitude ratios. On the critical isotherm the magnetic susceptibility 
behaves as 

X = CyH'-'/'; (2.14) 
the spontaneous magnetization vanishes as 

M = B{-tf. (2.15) 

One can then define the following universal ratio: 

Rc = aA+C+/B^, (2.16) 

which corresponds to the relation between exponents 

a + 2/3 + 7 = 2 . 

Indeed using this relation we verify that Rc is proportional to F{0, t)M~'^x which 
is normalization independent. Aharony and Hohenberg define the following uni- 
versal combination 

R^ = B^-^C+/{C^5Y , 

which corresponds to the relation ( |2.11|) . It is related to the quantity Qi defined 
by Fisher and Tarko, R^ = Qi^ . 

3 The e-expansion 

Let us first recall the results concerning exponents and the equation of state 
which have been obtained within the framework of the e = 4 — d expansion. 

3.1 Critical exponents 

Although the RG functions of the (p'^ theory and therefore the critical exponents 
are known up to five- loop order we give the expansions only up to the order 
needed in this article, i.e. to three loops, referring to the literature for details. 
The zero g*{e) of the /3-function is g*{£) = 167r^£/3 -|- 0{e'^). The values of the 
critical exponents 7, (3 and d = 1 + j / (3 are 

2/3 = l-i.+ ^.^+(^-Ac(3)).3 + +o(e^) 

in which C(s) is the Riemann ^-function. (At orders four and five C(5) and ^(7) 
successively appear. It has also been shown that starting at six loop order new 
non-(^ like numbers will appear |^2|.) Other exponents can be obtained from 
scaling relations. 
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3.2 The scaling equation of state 

The £-expansion of the scahng equation of state (eq. (|2.3|)) has been determined 
up to order £^ for the general 0{N) model, and order for = 1. Adjusting 
the normalizations of H and M to simplify the analytic expressions one finds at 
order 

fix) = l + x + efiix)+ e^f2ix) + e^fsix) + O {e^) , (3.1) 



with: 



with 



/i(x) = i(x + 3)L 

f2ix) = ^{x + 9)L2 + §^{x + 3)L 

L3 (27 + x) (675 + 246x + 25x2) 

■'^^^^ ~ 1296 ^ 1944 (3 + x) 

L (1617 + 701x - 1296xC(3)) x^ (-1 - 2A + 4C(3)) 
^ 17496 ^ 108 (3 + x) 

L = log(x + 3), 



while A is a constant 

A = \tjj'{l/3) - f TT^ = 1.17195361934 .... (3.2) 



The expression ( p.l|) is not valid for x large, i.e. for small magnetization M. In 
this regime the magnetic field H has a regular expansion in odd powers of M, 
i.e. in the variable z oc x~^ . Substituting in eq. ( ^Tj) x = xoz~^^^ (the constant 
xo takes care of the normalization of z) and expanding in e one finds at order 
for the function ( |2.6| ) 

F{z) = Fo{z)+eFi{z)+e^F2{z)+e^Fs{z), 

with 

Fo = z + ^z^ (3.3a) 
Fi = ^{-z^ + L{2z + z^)) 

F2 = {-50z^ + L{100z - Az^) + L^{18z + 27z^)) (3.36) 

^3 = 69^84 aT^^ \l\108z + 540z^ + 243z^) + L\1800z + U94z^ + 621z^) 
[2 + z ) L 

+ L(-1622;2^ + z^{1504: - 5184C(3)) + z(5608 - 10368C(3))) 

+z^ (-997 - 648A + 3888C(3)) + z^ (-2804 + 5184C(3))] (3.3c) 

and 

L = log(l + 2V2)- 
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While within the framework of the formal e-expansion one can easily pass 
from one expansion to the other one, still a matching problem arises if one wants 
to apply the ^-expansion for d = 3, i.e. e = 1. One is thus naturally led to 
look for a uniform representation of the equation of state valid in both limits. 
Josephson-Schofield parametric representation has this property. 

4 Parametric representation of the equation of state 

We parametrize M and t in terms of two variables R and 9, setting : 

M = moR^e , 

t = R{i-e'^), (4.1) 
H = hQRi^^h{e) , 

where ho, mo are two normalization constants ||3^]. This parametrization also 
corresponds in terms of the scaling variables x of eq. ( |3. 1| ) or z from eq. ( |2.4| ) to 
set 

z = p9/{l-9^f, e>0, (4.2) 
X = xop-'/^ (1 - 9^) 9-'/^, (4.3) 

where p is some other positive constant. 

Then the function h{9) is an odd function of 9 regular near 9 = 1, which is x 
small, and near ^ = which is x large. We choose ho such that 

h{9) = 9 + 0{9^) . 

The function h{9) vanishes for 9 = 9o which corresponds to the coexistence 
curve H = 0,T < Tc- From 9o we obtain the universal rescaled spontaneous 
magnetization \zo\ 

\zo\ = p9o/{9l-ir. (4.4) 

Note that the mapping ( f4.2|) is not invertible for values of 9 such that z'(9) = 0. 
One verifies that the derivative vanishes for 9^ = 1/(1 — 2/3) ~ 2.86. We shall 
see that this value is reasonably larger than the largest value of 9^ = 9q we will 
have to consider. 

Finally it is useful for later purpose to write more explicitly the relation be- 
tween the function F{z) of eq. ( p.6|) and the function h{9): 

h{9) = p-^ {1 - 9^Y^ F {z{9)) . (4.5) 
Expanding both functions 

F(z)=^+i;^3^5^F2;+i^2Z+i^ (4.6) 

1=2 

h{9)/9 = l + J2h2i+i9^', (4.7) 
1=1 
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we find the relations 



hs = y-^ (4.8a) 

/i5 = h(7 - 1) + i(2/3 - 7)P' + Fsp^ (4.86) 
hr = ^7(7 - 1)(7 - 2) + 1^(2/3 - 7)(2/3 - 7 + l)p2 

+ (4/? - ^)Fry + Fjp^ (4.8c) 



From the parametric representation of the equation of state it is also possible to 
derive a representation for the singular part of the free energy per unit volume. 
Setting: 

J-sg.(M,t) = rsg.(M,t)/vol. = homoR^-''g{e), (4.9) 
one finds for g{9) a differential equation: 

h{9) (1 - 6*2 + 2/36*2) = 2(2 - a)9g{9) + (l - 9^) g'{9). (4.10) 

The integration constant is fixed by requiring the regularity of 5f(^) at ^ = 1. In 
the same way the inverse magnetic susceptibility is given by: 

X"' = {ho/mo)R''g2{9), (4.11) 



with: 

g2{9) {1-9^ + 2/36*2) = 2p69h{9) + (l - 9^) h'{9). (4.12) 

Note that g2 then in general has a pole dX 9"^ = — 2/3), the point at which 
the mapping ( [4.2| ) is not invertible. 

The functions g{9) and g2{9) will be used to calculate several universal ratios 
of amplitudes. 



4-1 Amplitude ratios 

The amplitude ratios defined in section ^ are related to the scaling equation of 
state. They can be calculated from the functions h{9), g{9) and (72(6*) appearing 
in the parametric representation, and thus ultimately from the function h{9) 
alone (eqs. ( [4.10| , |4.1^ )). One verifies that indeed they do not depend on the 



variable R and the constants mo, Hq appearing in ( fl.l|) . 

The magnetic susceptibility. The magnetic susceptibility in zero field can be 
calculated from the function g2{9) defined by equation ( }4.12| ). One obtains (equa- 
tion (|T2D): 

^+_..2 , h'{9o){l-9l) 

C^"^^° ^> h'{0){l-9l + 2(39^,y ^^■^'^> 
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The specific heat. This ratio is directly related to the function g{9) defined by 
equation ( [4.1C1|) : 

Other universal ratios. Similarly the ratios Rc and can be derived from 
the functions g{9) and h{9) 



h'{oM-i] 



= _ «)(2 - a f' ZLn2 ^ (4-15) 



= J:\/2 (4-16) 



4.2 Parametric representation in the e-expansion 

Up to order e"^ the constant mo (or p) can be chosen in such a way that the 
function h{9) reduces to: 

h{e) = e {1 - le^) + o {e^) . (4.17) 

The simple model in which h{9) is approximated by a cubic odd function of 9 
is called the linear parametric model. At order e'^ the linear parametric model 
is exact, but at order the introduction of a term proportional to 9^ becomes 
necessary [^,23]. One finds: 

h{9) = 9{1 + hs9^ + h^9^) + O {e^) , (4.19) 

with 

^3 = -^(l + f^), = 1^ (C(3) - iA - i) , (4.20) 



where A is the constant ( p.2|) . 

The function h[9) vanishes on the coexistence curve for 9 = 9q: 

^o'-^(l-f^)+0(.^)- (4.21) 

Note that /13 and thus 9q are determined only up to order . It follows 

p2 = 6(7 + /la) = 2 (1 + i£ + ^e^) = 3.13 ± 0.13, 

because /la is determined only up to order , and 

W\=p9o{9l-l)-f' 

= V3 X 2/^ [1 + i£ + + (^A - ^C(3) + Sk) e'] ~ 2.87 ± 0.06 



12 



(in this case we summed by a Pade [1,2]). 

Remark. Even in the more general 0{N) case, the parametric representation 
automaticaUy satisfies the different requirements about the regularity proper- 
ties of the equation of state and leads to uniform approximations. However for 

> 1 the function h{9) still has a singularity on the coexistence curve, due to 
the presence of Goldstone modes in the ordered phase and has therefore a more 
complicated form. The nature of this singularity can be obtained from the study 
of the non-linear cr-model. It is not clear whether a simple polynomial approx- 
imation would be useful. For = 1 instead, one expects at most an essential 
singularity on the coexistence curve, due to barrier penetration, which is much 
weaker and non-perturbative in the small e- or small (^-expansion. 

Amplitude ratios. The use of eqs. ( [4.14| - ^T6|) together with the expression of 
h{9) eqs. (^4.19| , [4.20|) evaluated at £ = 1 gives us the predictions for the amplitude 
ratios reported in table 5 (e-expansion (6)). Moreover from these equations the 
known e-expansion of various amplitude ratios, ||2^ , p3| , |25[| , can be easily obtained. 
The ratio C+/C~, related to the magnetic susceptibility in zero field, is 



C- 6/3-1 - > ^ ' •> ^''t^') 

= 2-'-'(* - 1) [1 + (C(3) + fA + i) £'] (4.226) 

The ratio /C~ can be expressed at order entirely in terms of critical ex- 
ponents. This form follows naturally from the parametric representation of the 
equation of state. The relative correction is of the order of only 3%. The 
three first expressions yield for e = 1 respectively (the exponents being replaced 
by the central values of the summed e-expansion) : 

C+ 

— = 4.688, 4.757, 4.863 . 

The spread gives for this short series an indication about the uncertainty about 
the result. 

The ratio A'^ /A~ at order £^ is given by 

^ = 2-^ [l + .-f(i-iA-C(3)).^]+0(.3). 
The ^-expansion of Rc is: 

Rc = y^'-'e [1 + + - |C(3) - fA) e^] + O {e') . (4.23) 
R^ is given by 

R^ = 3(^-3)/227+(i-5)/2 [1 + + - ^C(3)) + O {e') . (4.24) 
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5 The perturbative expansion at fixed dimension three 

Critical exponents and several other universal quantities have been estimated 
using the perturbative expansion at fixed dimension d < 4. Since the massless 
theory is then IR divergent, calculations have been performed within the frame- 
work of the 3D perturbative expansion renormalized at zero momentum for the 
massive cp'^ field theory. As in the case of the ^-expansion, it is necessary to first 
determine the IR stable zero g* of the function (3{g) which is given by a few 
terms of a divergent expansion. The obvious problem is that we have no longer 
a small parameter in which to expand and g* is a number of order 1. Already 
at this point a summation method is required. Note also that at any finite order 
the results for universal quantities become renormalization scheme dependent in 
contrast with the results of the e-expansion. 

On the other hand, because one-loop diagrams have, in 3 dimensions, a simple 
analytic expression, it has been possible to extend the calculation of RG functions 
in the A^-vector model up to six-loop order. The expansions are [T 



f3{g) = -g + f- 111^3 ^ 0.3510695978^^ - 0.3765268283^^ 

+ 0.49554751^^ - 0.749689^^ + O (f) , (5.1) 
-f-^{g) = 1 - 1^ + ±g^ - 0.0230696213^^ + 0.0198868203^^ 

- 0.02245952^^ + 0.0303679^^ (5.2) 
rj{g) = + 0.0009142223^^ + 0.0017962229^^ 

- 0.00065370^^ + 0.0013878^^ (5.3) 

with the normalization: 

~g = 3^/(167r). (5.4) 

It has thus been possible to determine critical exponents quite accurately (table 
1)- 

To determine the equation of state or universal ratios of amplitudes a new 
problem arises. In this framework it is more difficult to calculate physical quan- 
tities in the ordered phase because the theory is parametrized in terms of the 
disordered phase correlation length ^ = which is singular at Tc (as well 

as correlation function normalization condition, eq.( |2.1| )). Let us consider the 
perturbative expansion of the scaling equation of state ( |2.6| ) . For example at 
one-loop order for d = 2> the function F{z,g) is given by: 

F(z, g) = z + lz^- ^gz [(1 + zV2)'^' - 1 - ^V4] (5.5) 

= z + \z^ + ^gz' - -^gz' + 0(z9), 

where the subtractions, due to the mass and coupling normalizations, are deter- 
mined by the conditions (|2.1j). This expression is adequate for the description 
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of the disordered phase, but aU terms in the loopwise expansion become singu- 
lar when t goes to zero for fixed magnetization, that corresponds to the hmit 
m — > 0, fixed for the effective potential, or 2; ^ oo for reduced variables (see 
eq. (|2.4|) ). In this regime we know from eg . (12 .31 ) (or from RG considerations) 
that the equation of state is 

H{M, t = 0)(xM^ F{z) oc (5.6) 

exactly. 

In the case of of the e-expansion the very essence of the method is that one 
is doing perturbative expansions for the theory at the critical point: it follows 
that the scaling relations (and thus the limit eq. (p.6|)) are exactly satisfied order 
by order. Moreover the change to the variable x oc z~^^^ (more appropriate for 
the regime t — 0) gives an expression for f{x) oc F{x~^)x^^ that is explicitly 
regular in a; = (Griffith's analyticity at the critical point): the singular powers 
of logx induced by the change of variables cancel non trivially at each order, 
leaving only regular corrections. 

The situation changes when one deals with the perturbation theory aX d = 
3 dimensions: for g generic the system is no more at the critical point, and 
consequently scaling properties are not satisfied order by order in g. In particular 
the change to the Widom function f{x) will introduce the singular terms in logx 
that violates Griffith's analyticity. An analogous problem arises if one first sums 
the series dX g = g* before changing to the variable x. In this case the singular 
contributions (in the form of powers of x) do not cancel, as a result of unavoidable 
numerical summation errors. 

Several approaches can be used to deal with the problem of reaching the or- 
dered phase. In [|2^ a method to calculate amplitude ratios is proposed which 



has been also used in | 

If we are concerned only with amplitude ratios another strategy is available 
that in some sense bypasses the problem. Near the critical point physical quan- 
tities have simple power law singularities in t. To reach the coexistence curve, 
i.e. t < 0, -ff = 0, it is possible to proceed by analytic continuation in the complex 
t-plane. From equations ( |2.8|J2.10|) 



and the knowledge that at M fixed H{M, t) is regular at Tc or t = we conclude 
that t < corresponds to z complex 

t = Itle'"", ^ z = Izle-'""^ . (5.7) 

The scaling variable H(—t)~^^ is then given by: 

H{-t)-^^ = e'-^l^^ F{z) = \F{z) I . (5.8) 
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Finally the spontaneous magnetization is given in terms of the complex zero zq 
ofF{z). 

It is in particular possible to evaluate ratios of amplitudes of singularities above 
and below Tc. we calculate the complex zero zo{g) of F{z,g) and substitute 
it in other quantities. The result is complex but its absolute value converges 
towards the correct result. For example the ratio of amplitudes for the magnetic 
susceptibility (eq. (|2.10| ) ) is given by: 

C+/C- = e— ^ F' {zo (g) , g) /F' (0, ^7) = \F' {zo (g) , ^7)! • (5.9) 

We thus get a series expansion in g for /C~ which can be summed at g = g* 
with techniques described in section ^ However this method does not allow to 
calculate for t small and thus is not well suited to determine the full equation of 
state. Moreover if we want to sum perturbation series with the method recalled 
in section we will face the problem that the large order behaviour depends on 
the value of the variable z itself. 

Encouraged by the results obtained within the ^-expansion scheme, we develop 
a different, more powerful strategy, based on the parametric representation. 

The parametric representation. Order dependent mapping (ODM). The prob- 
lem that we face is the following: to reach the ordered region t < for the 
(summed) equation of state function F{z, g*), we must cross the point 2; = 00 for 
which the exact behaviour is dictated by eq. ( p.6|) . The idea that the e-expansion 
suggests is to introduce an new field variable 9 and an auxiliary function h{9) 
defined as in ([4.2| , [4.5|) : in this way the exact function h{9) will be regular near 
9 = 1 (i.e. z = 00) and up to the coexistence curve. However, the approximate 
h{9) that we obtain by summing perturbation theory at fixed dimension, will still 
not be regular. The singular terms generated by the the mapping eq. ( [4.2| ) at 
9 = 1 will not cancel exactly due to summation errors. The last step we propose 
is to Taylor expand the approximate expression of h{9) around 9 = and to 
truncate the expansion, enforcing in this way regularity. The next question then 
is to which order should we expand? Since the coefficients of the 9 expansion 
are in one to one correspondence with the coefficients of small z expansion of 
the function F{z,g*) the maximal power of 9 in h{9), should be equal to the 
maximal power of z which can be reasonably well summed. Indeed although 
the small z expansion of F{z) at each finite loop order in g contains an infinite 
number of terms, the determination of the coefficients of the higher powers of z 
is increasingly difficult. The reasons are twofold: 

(i) The number of terms of the series in g required to get an accurate estimate 
of Fi increases with / (see section P) . 

(ii) At any finite order in g the function F{z) has spurious singularities in the 
complex z plane (see e.g. eq. (|5.5|) , z'^ = —2) that dominate the behaviour of the 



coefficients Fi for / large. 

In view of these difficulties we have to ensure the fastest possible convergence 
of the small 9 expansion. For this purpose we use the freedom in the choice 
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of the arbitrary parameter p in eq.( [4.2|) : we determine it to minimize the last 
term in the truncated smaU 9 expansion, i.e. increasing the importance of smaU 
powers of theta which are more accurately determined. This is nothing but the 
application to this particular example of the series summation method based on 
ODM PI. 



This strategy applied to the available data, leads at leading order for h{9) to 
a polynomial of degree 5, whose coefficients are given by the relations 



h{e) = e[i + /i3(p)^' + h^{p)e% (s.io) 



For the range of admissible values for F5 the coefficient of 6^ given by eq. (|4.8fc| ) 
has no real zero in p. It has a minimum instead 

p' = P5 = Y^(7-2/3). (5.11) 

Substituting this value of p into expression ( p. 101 ) we obtain the first approxi- 
mation for h{6). At next order we look for a minimum p-j of |/i7(p)|. We find a 
polynomial either of degree 5 in ^, when /ly has a real zero, or of degree 7 when 
it has only a minimum. 

We have not explored beyond h^^p) because already Fg is too poorly deter- 
mined. 

Finally we want to point out that, to compute the amplitude ratios and thus 
to reach the ordered phase quantities we implement in the present framework the 
idea of analytic continuation outlined in the previous subsection. In particular 
the ordered phase at zero magnetic field will correspond to the choice 6 = 6q 
where 9q will be the zero of h{6) closest to the origin. 

6 The equation of state: series summation 

Our first task is thus to determine the coefficients F21+1 as accurately as possible. 
To sum the perturbation series the Borel-Leroy transformation has been 
used, combined with a conformal mapping |^] (a simplified version of method 



used in for critical exponents). Let S{g) be the function whose series has to 
be summed. We then transform the series: 

S{g)=Y,Skg\ (6.1) 

fc=0 

into: 

/"OO 

Sig) = TBk{b) t^e-'u\gt)dt, (6.2) 

with: 



VI + as - 1 
Vl + as -f- 1 
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The coefficients Bk are calculated by expanding in powers of g the r.li.s. of 
equation ( |6?^ ) and identifying witli expansion ( |6.1j ) . (For motivations and details 
see e.g. [|1|). The constant a (corresponding to the scaled variable g = -^g of 
@) has been determined by the large order behaviour analysis 

a{d = 3) = 0.147774232 , (6.4) 

for the perturbative expansion in d = 3 dimensions. The parameter b is adjusted 
empirically to improve the convergence of the transformed series: for example 
one looks for the intersection between the results at two consecutive orders in k 
(or the minimum of the difference) , or a point of least sensitivity. Moreover the 
value of b has to stay in a reasonable range around the value predicted by the 
large order behaviour. It is also to be expected that the summation method will 
be efficient if the coefficients Sk are already approaching the asymptotic large 
order regime. This is what we are going to test first. The expectations are that 
for the coefficients F21+1 the situation will deteriorate with increasing /: indeed 
the large order behaviour estimate has increasingly large corrections. Moreover 
for / large F21+1 becomes dominated by the perturbative singularity of F{z) at 
z'^ = —2 while the summed function has singularities at different locations. 



Large order behaviour analysis The series are taken from p6| , p7[|28[] . The 
coefficients of the perturbative expansion of 

Fi ^ i^7i-(^+i)/2m('+^)/2-3r('+i)(p, = 0,m,^7) = J^F^./, 

k 

have the large order behaviour: 



Fi, ~ (-A)^-i+(m)/2r(fc + / + 3/2) 

g-974/(327r) /4/2y/2 

Cl = 



TC 



(27r)3/2 1^ /4 



3 T 



-1/2 



The constant I/^ = (a = a{d = 3) is given by eq. (|6^)), is related to the 
action of the instanton which describes the instability of the (p"^ for negative 
coupling, while Ii = 31.691522, Iq — 659.868352, (related to integrals of power 
1, 6 of the scaled instanton solution) and Di = 10.544 (related to the functional 
determinant) have been computed in 



Table 3 shows that, as anticipated, the asymptotic regime sets in later when / 
increases. We thus expect that the efficiency of the summation will correspond- 
ingly decrease, and indeed this is what happens. 

Summation. Following an idea introduced in ^ for the summation of the 
^-expansion we have in addition made a homographic transformation on the 
coupling constant g to displace some singularities in the complex ^r-plane: 

g = g'/il+pg'/gn- (6-5) 
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Table 3 

Large Order Behaviour Analysis : Fik/ Fj^^ . 



I 


k = 2 


fc = 3 


fc = 4 


= 5 


5 


1.9646 


1.5894 


1.4226 


1.3498 


7 


1.9680 


1.9965 


1.9603 


1.9204 


9 


1.4443 


1.9879 


2.229 


2.333 



We have looked for a value of the parameters p and h for which the results 
were speciaUy insensitive to the order k: in practice we minimized the absolute 
differences of results corresponding to three successive orders. Finally to verify 
that our method gives results consistent with those of the previous analysis we 
have applied it to the RG /3-function and exponents. Figure 1 gives the last four 
orders for the /9-function as a function of the parameter h and at the minimal 
value p = 0.196. One observes that the curves flatten with increasing order, as 
expected. The central value slightly differs from the result given in 

g* = 23.70 ± 0.05 , with u = (3'{g*) = 0.79 ± 0.01 , 

and the apparent error is smaller. Note that a variation of g* of 0.05 yields a 
variation of 0.04 of the /^-function. 

Figure 2 displays the results for the exponent 7 for the four last orders, for 
g* = 23.70, and for the minimal value p = .182. One obtains 

7 = 1.2405 ±0.0012 . 

With the same method one finds for the exponent /3 

(3 = 0.3250 ±0.0015 . 

One notices the prefect agreement with the result of table 1 , although the central 
value for g* has changed. The reason is that the central values given in for 
the exponents largely rely on the pseudo-£ expansion, which avoids the explicit 
determination of g* . 

7 Numerical results: Equation of state 

We first calculate the first coefficients of the small z expansion of the equation of 
state function F{z) (or equivalently of the reduced effective potential, eq. (pT6|), 
summed at the I.R. fixed point g*). In figure 3 we display the behaviour of -F5 
in terms of b at the minimal value of the parameter p = .182 (see section |^ for 
definition of parameters). 
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Fig. 1 The RG /9-function plotted vs. the parameter h for successive orders k. 



Exponent Gamma 
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1.242 
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k= 6 

— k=5 
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Fig. 2 The exponent 7 plotted vs. the parameter h for successive orders k. 

Table 4 contains our results together with other published estimates of the 
coefficients of the small z expansion of Fi^z). For what concerns the ^-expansion, 
we have taken the e-expansion of h{Q\ eqs. ( ^4.iy| . |ir^DD setting e = \. We have 
then calculated the coefficients -F2Z+i5 using the numerical values of the more 
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Coefficient F 5 
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Fig. 3 The coefficient F5 plotted vs. the parameter b for successive orders fc. 



accurately determined exponents 7, /3 (second line of table 1). This procedure 
has led to much stabler results for F21+1 than a direct summation of the e- 
expansion. The quoted errors are nevertheless only indicative (and subjective) 
because h{9) is known only to order £^ for /i3 and for h^. Note that although 
g* is known only to order £^ EOl 



g%e) = ir(rf/2)(47r)^/^£ (l + §5) + 0(8^), 



evaluation of the expression at £ = 1 gives nevertheless a number of the right 
order of magnitude g* = 28. 

Parametric representation. We then determine by the ODM method the coef- 
ficient p and the function h{9), as explained in section ^ We obtain successive 
approximations in the form of polynomials of increasing degree for h{9). Note 
that we here have a simple test of the relevance of the ODM method. Indeed, 
once h{9) is determined, assuming the values of the critical exponents 7 and 
we can recover a function F{z) which has an expansion to all orders in z. 
As a result we obtain a prediction for the coefficients F21+1 which have not yet 
been taken into account to determine h{9). The relative difference between the 
predicted value and the calculated one gives an idea about the accuracy of the 
ODM method. Indeed from the values F5 = 0.01712,7 = 1.2405, /3 = 0.3250, we 
obtain 



F7 = 4.94 X 10"^, Fg = -3.2 x 10-^ Fu = 8.5 x 10" 
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Table 4 
Equation of state. 



9 



F7 



Fg 



e— exp. 
d = 3 
d = 3 |1| 

HT m 
HT m 
HT m 

MC |4| 
MC [Hi 
ERG [|6| 
ERG |7| 



28. 

23.70 ±0.05 

23.73 
23.72 ± 1.49 
24.45 ±0.15 
23.69 
23.3 ±0.5 
24.5 ± .2 

28.9 
20.7 ± 0.2 



0.0176 ±0.0004 
0.01712 ±0.00006 

.01727 
0.0205 ±0.0052 
017974 ± .00015 

.0166 
0.0227 ±0.0026 
0.027 ±0.002 

0.016 
0.0173 ±0.0001 



(4.5 ±0.3) X 10"'' 
(4.96 ±0.49) X 10"'' 
.0010 



.00055 



0.00236 ± .00040 

4.3 X 10"'' 
(5.0 ±0.2) X 10"" 



-(3.2 ±0.2) X 10"^ 
-(6.2 ±4.0) X 10"^ 



.00001 



-(4. ±2.) 



10 



-5 



We note that the value for F7 is very close to the central value we find by direct 
series summation, while the value for Fg is within the error bars. This result gives 
us confidence in our method. It also shows that the value of Fg obtained by direct 
summation contains very little new information, it provides only a consistency 
check. Therefore the simplest representation of the equation of state, consistent 
with all data, is given by 

h{e)/e = 1 - 0.76147(36) + 7.74(11) x 10"^ ^^ (7.1) 

(errors on the last digits in parentheses) that is obtained from = 2.8743 (fixed 
according eq. ( |5.11|) ). This expression of h{9) has a zero at 

el = 1.331, (7.2) 

to which corresponds the value of the complex root zq of F{z), \zo \ = 2.801 (the 
phase, given by eq. (|5.7|) , is —in (3). The coefficient of 9^ in eq. ( |7.1|) is smaller 
than 10""^. Note that for the largest value of 9^ which corresponds to 9q, the 
9"^ term is still a small correction. Finally note that the corresponding values 
for the ^-expansion are hs = —0.72, = 0.013. These values are reasonably 
consistent between them, since a small change in the value of p with a correlated 
change in induces a very small change in physical quantities. 

The Widom scaling function /(x), eq. (pT3D, (with /(—I) = and /(O) = 1) 
can easily be derived by (numerically) solving the following system: 

fix)=9-'hi9)/h{l) 
X = ■ 



91J\9 
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Our results are displayed in fig. 4, where they are compared with £^ (Linear 
Parametric Model) and predictions, High Temperature series ("quintic" fit 
amplitude ratios results of |^ and Pade approximants of bcc lattice results 
of 1^^) as well with Exact Renormalization Group (a fit reported in ||5^ that 
approximates numerical data up to 1%). Note that the rapid increase of relative 
errors for x < is not very significative since the absolute values are small. 
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Fig. 4 Widom scaling function fix): relative errors with respect 3D QFT plotted 
vs. X. 

One notices the general agreement between different predictions within 6% up 
to a; = 10. In particular the general agreement between our predictions and the e- 
expansion is striking. The main disagreement with other predictions comes from 
the region x oo, i.e. from the small magnetization region, where our predictions 
should be specially reliable. Finally let us notice that the results are based 

some obsolete estimates of critical amplitudes (7 = 5/4, /3 = 5/16). 

8 Amplitude ratios. Conclusion 

We finally use the expressions ( [1.14j - pT6|) with the critical exponents determined 
by earlier work (from longer series) to calculate various amplitude ratios. As a 
comparison we apply the same procedure to the e-expansion at order (using 
the values of exponents determined from the ^-expansion at order e^). It is 
necessary to again stress that as in table 4 the errors we then quote are largely 
tentative since the order is the first order where some reliable results may be 
expected. 
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Table 5 contains a comparison of amplitude ratios as obtained from RG, lat- 
tice calculations and experiments on binary mixtures, liquid- vapour, uniaxial 
magnetic systems. 



Table 5 
Amplitude ratios. 





A+/A- 


C+/C- 


Rc 


Rx 


£ — exp. (a) 


0.524 ±0.010 


4.9 




1.67 


£ — exp. (6) 


0.547 ±0.021 


4.70 ±0.10 


0.0585 ±0.0020 


1.649 ±0.021 


d = 3 fixed (a) 


0.541 ±0.014 


4.77 ±0.30 


0.0594 ±0.001 


1.7 


d = 3 fixed (6) 


0.536 ±0.019 


4.82 ±0.10 


0.0576 ±0.0020 


1.674 ±0.019 


HT series 


0.523 ±0.009 


4.95 ±0.15 


0.0581 ±0.0010 


1.75 


bin. mix. 


0.56 ±0.02 


4.3 ±0.3 


0.050 ±0.015 


1.75 ±0.30 


liqu. — vap. 


0.48-0.53 


4.8-5.2 


0.047 ±0.010 


1.69 ±0.14 


magn. syst. 


0.49-0.54 


4.9 ±0.5 







Results for e-expansion (a) are taken from ||2^ and (direct Fade summation 
of each corresponding series), while (6) are our results, obtained by first summing 
h{6) and then computing ratios by use of eqs. The results (i = 3 

fixed dimension (a) are taken from [|2^ and refers to direct summation up to 
0{g^) while d — 3 (6) is the present work. High Temperature results are taken 
from {R-^ from Experimental data are extracted from |5^], to which 

we refer for an updated and wider list of results (and references). We note the 
general consistency of the results obtained by different methods. In addition 



in table 7 we compare our predictions with the results of concerning the 
following the universal ratios: 



Rq 



R? 



91 

C7 



where the are defined by d^x/dH^ = C^_^2W other quantities 

have been defined in section |^. After some algebra one finds the relation 

Ro = —(1 - 12^5)-^ 



that has thus been used in Table 4. 



Conclusion. Working within the framework of renormalized quantum field 
theory and renormalization group, we have shown that the presently available five 
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Table 6 
Other amplitude ratios. 







Rs 




HT series |3S] 

d=3 fixed 
e-expansion 


0.1275 ±0.0003 
0.12586 ±0.00012 
0.1268 ±0.0008 


6.4 ±0.2 

6.10 ±0.06 

6.11 ±0.10 


-9.0 ±0.3 

-9.2 ±0.6 
-8.3 ± 1.0 



loop series allow, after proper summation, to determine with reasonable accuracy 
the complete scaling equation of state for 3D Ising-like systems. The parametric 
representation of the equation of state plays a central role in our analysis. As a 
consequence new estimates of some amplitude ratios have been obtained. Clearly 
a similar strategy could be applied to other quantities in a magnetic field, in the 
scaling region. We want to stress that an extension of the e-expansion of the 
equation of state for = 1 to order or even better e^, that does not seem 
an impossible task, would significantly improve the e-expansion estimates and 
would therefore be extremely useful. Moreover the present approach could be 
extended to systems in the universality class of the {(p'^)s field theory for higher 
N, provided expansions of the renormalized effective potential at high enough 
order are computed. 
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APPENDICES 



Al Perturbative expansion of the eflfective potential 



Table 7 
Coefficients Fuk- 



h 


I 


A; = 


k = 1 


= 2 


1 





-l/(127r) 






2 





1/(12871^) 






2 


1 




-l/(967r2) 




3 





(5 + 81og|)/(6 1447r3) 


l/(7687r3) 




3 


1 


A;3/(69127r^) 






3 


2 


A;4/(20 7367r^) 






4 





0.449205291 x 10"^ 


0.122593698 x lO"'^ 




4 


1 


-.115567060 X 10"^ 


-l/(184327r'^) 




4 


2 


0.100772569 x lO"'^ 






4 


3 


-.300586776 x 10"*^ 






5 





-.451638912 x 10~^ 


(3-81ogf)/(11796487r^) 


-l/(29 49 127r^) 


5 


1 


0.806878091 x 10"^ 


^3/(132710471-^) 




5 


2 


-.100671980 X 10"*^ 


^4/(132710471-^) 




5 


3 


0.647972392 x lO"'^ 






5 


4 


-.165599191 X lO"'^ 







Our starting point is the renormalized five loop eiTective potential in Minimal 
Subtraction (MS) scheme, that can be written in the form: 

TTl ^ A 

6=1 1=1 k=0 ^ ^ \ / ^ ^ 

where (A4) are the renormalized mass (coupling), and we defined = 
rus'^ + ^0^, ki = y/e/3. The constant F^k are reported in Table 7. 
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Analytical results up to three loops have been obtained by [^Sf, and can be 
parametrized by the two constants 

^3 = -| + ^ - y (log(4/3))' - 27Li2(l/4) 
= 9.36271548728022813982 ■ ■ ■ 



^ + ^(log(4/3))' + 541og(4/3) + yLi2(l/4) 

54 /"^ dx / 3 , 3 + x x^ , 4 x , 3 + x 

log - + log — 5" log — h — — log — - — 



72 Jo V 4 ^2 + x 4-x2 ^2 + x 2 + x 

= -6.43307044049269064141 • • • 

These expressions involve the dilogarithm or Spence function 

We have obtained the analytic expressions of some four and five loop coefficients 
by imposing the MS RG equation on the effective potential, that in our notations 
reads: ^ 

+ WT^ — TT-) ^MS = 0- (^1-2) 

oil 12{4n)^ nis drus J ^^^^ 



We want to emphasize that eq. ([A1.2| ) above holds to all orders in this scheme 

because mass divergences arise only at two loop order (A2 = + i2(47r)2 

and no additional finite renormalization is needed in the MS scheme. In prac- 
tice ( |A1.2|) should be satisfied at b loops up to order O^X^') (if we take (p ~ 
0{1/ the term enforces relations between the coefficients of log, log^ 
at order b and terms at order 6 — 2. The constant -^4,0,1 cannot be fixed be- 
cause the corresponding relation implied by ( |A1.2|) involves terms without field 
dependence (i.e. "cosmological terms") that are not considered in ( [Al.l| ). 



Numerical results up to five loops for Fi,ik have been first published in Table II 
of pH] (derived from Table I that contains the values of single diagrams, see also 



Some misprints in Table I, of eq. (3.2) and Table II of |26] are reported in 



where Table II is rederived from the corrected values of Table I, a corrected 



eq. (3.2) and eq. (B1-B3) of [|2^. We checked once more Table 1 of [|2^, and we 
obtained the same values, provided that — (— A^)^ of eq. (3.2) of is replaced 
by {—Ad}^ (a misprint is present in the prescriptions of [^), provided that 
the formulae eqs. (B1-B3) of PH] are intended to hold for the dimensionless Ti, 



functions (i.e. divided by a factor Qq compared to definition eq. (3.2)), and that 

d 

dro 



Xq is considered as dependent from tq when taking derivatives ^J- as in eq. (B3). 
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For our purposes we are interested in the physical mass scheme, that is fixed 
by normahzation conditions ( |2.1| ). The required finite renormahzation can be 
summarized by the foUowing relation: 

Vm((/), m, g) = VMsi\^(p, Zmm, ZgUig, Jj) (Al.3) 

that defines the effective potential in the new scheme in terms of the new renor- 
malized mass m and dimensionless coupling g. 

The renormahzation constant Z is fixed by (^T^) and we rewrite here for com- 
pleteness. It is obtained from the derivative of the two-point function: 

/ p=0 \ / p=0 

The last equality is due to the absence of field renormahzation between bare and 
MS quantities. Thus the constant Z is the same that one would obtain in passing 
from bare quantities to mass scheme and can be obtained from the known RG 
functions of the latter scheme: 



JO M 

The other constant can then be fixed by 

— ^(0, m,5f) = m 



-^^(0,m,^) = mg. 

Due to the fact that the potential in the mass scheme will not depend on //, 
one can simplify tedious calculations by choosing Ji oc m^. 
The reduced potential V is then obtained from: 

Tfl J 

y{z,g) = — VM{z^Jm/g,m,g). 
9 

We report here the final (numerical) results for the coefficients Fi of the Taylor 
expansion in z of F{z^ g), the derivative with respect to 2; of V (see (|2.6| )). 

^ ^ Q_ _ g" 9^ (-274 + 18A;3-24A;4 + 271og(|)) 
^ 2567r 20487r2 4423687r3 

- 2.25018021 X 10" V + 1-977252 x 10"V 

g 65^^ g^ (646 - AQk^ + 80^4 - 45 log(|)) 

^ ~ ~ 10247r ^ 1474567r2 ^ 17694727r3 
+ 4.17878698 x 10"^/ - 4.512966 x lQ~^g^ 

5g 25g^ g^ (-69099 + 4200A;3 - 11200/c4 + 3780 log(f )) 

^ ~ ^16384^ ~ 983047r2 ^ 339738624^^3 

- 5.21924201 X lO'^^;'^ + 6.986358 x lO'^^^^ 
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